Phase Diagram of Colloidal Hard Superballs: from Cubes via Spheres to Octahedra 
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The phase diagram of colloidal hard superballs, of which the shape interpolates between cubes and 
octahedra via spheres, is determined by free-energy calculations in Monte Carlo simulations. We 
discover not only a stable face-centered cubic (fee) plastic crystal phase for near-spherical particles, 
but also a stable body-centered cubic (bcc) plastic crystal close to the octahedron shape. Moreover, 
coexistence of these two plastic crystals is observed with a substantial density gap. The plastic fee 
and bcc crystals are, however, both unstable in the cube and octahedron limit, suggesting that the 
rounded corners of superballs play an important role in stablizing the rotator phases. In addition, 
we observe a two-step melting phenomenon for hard octahedra, in which the Minkowski crystal 
melts into a metastable bcc plastic crystal before melting into the fluid phase. 
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Recent breakthroughs in particle synthesis have re- 
sulted in a spectacular variety of anisotropic nanopar- 
ticles such as cubes, octapods, tetrapods, octahedra, ice- 
cones, etc. [ij. A natural starting point to study the self- 
assembled structures of these colloidal building blocks is 
to view them as hard particles [1]. Not only can these 
hard-particle models be used to predict properties of suit- 
able experimental systems, but such models also provide 
a stepping stone towards systems where soft interactions 
play a role Moreover, the analysis of hard particles 

is of fundamental relevance and raises problems that in- 
fluence fields as diverse as (soft) condensed matter 0,0- 
0| , mathematics [1, 0] , and computer science [1] . In this 
light the concurrent boom in simulation studies of hard 
anisotropic particles is not surprising 

The best-known hard-particle system consists of hard 
spheres, which freeze into close-packed hexagonal (cph) 
crystal structures 8], of which the ABC-stacked cph 
crystal, better known as the face-centered cubic (fee) 
crystal phase, is thermodynamically stable [17i]. Hard 
anisotropic particles can form liquid-crystalline equilib- 
rium states if they are sufficiently rod- or disclike [l2l.[l6|. 
but particles with shapes that are close-to-spherical tend 
to order into plastic crystal phases, also known as rota- 



tor phases |14l - lla |. In fact, simple guidelines were re- 
cently proposed to predict the plastic- and liquid-crystal 
formation only on the basis of rotational symmetry and 
shape anisotropy of hard polyhedra . In this Letter we 
will take a different approach, based on free-energy cal- 
culations, and address the question whether and to what 
extent rounding the corners and faces of polyhedral par- 
ticles affects the phase behavior. Such curvature effects 
are of direct relevance to experimental systems, in which 
sterically and charged stabilised particles can often not 
be considered as perfectly flat-faced and sharp-edged [3] ■ 
For instance, recent experiments on nanocube assemblies 
show a continuous phase transformation between simple 



cubic and rhombohedral phases by increasing the ligand 
thickness and hence the particle sphericity l3| . 

In this Letter, we study a system of colloidal hard su- 
perballs in order to address these problems. A superball 
is defined by the inequality 



+ < 1, 



(1) 



where x, y and z are scaled Cartesian coordinates with 
q the deformation parameter. The shape of the super- 
ball interpolates smoothly between two Platonic solids, 
namely the octahedron [q = 0.5) and the cube [q — oo) 
via the sphere {q — 1). By determining the phase dia- 
gram of these superballs as a function of g, we discovered 
a thermodynamically stable body-centered cubic (bcc) 
plastic crystal phase for octahedron-like superballs. To 
the best of our knowledge no plastic crystals other than 
cph structures have so far been observed for hard parti- 
cles. Moreover, we find that bcc and fee plastic crystal 
phases are unstable for hard octahedra and hard cubes, 
respectively. Therefore, rounded faces and edges may 
play an important role in stabilizing rotator phases, while 
flat faces tend to stabilize crystals. 

Following Refs. [l^, we first calculate the close- 
packed structures for systems of hard superballs. We 
employ the algorithm of checking overlap described in 
Ref. [20|. For cube- like particles, it is found that at close 
packing there are so-called Co and Ci crystal phases 
in accordance with Ref. When we perform NPT 

Monte Carlo simulations with variable box shape to de- 
termine the equation of state (EOS) of the crystal phase, 
our simulation results show that both the Cq and the Ci 
crystals deform with decreasing density. The lattice vec- 
tors for Ci crystals are given by ei = 2^^^/^'' i + j 
,62 = 21-1/2? J ^ 21-1/29 k , 63 = 2(s + 2-^/^1) i - 2s j - 
2s k, where i, j and k are the unit vectors along the axes 
of the particle, s is the smallest positive root of the equa- 
tion (s + 2-1/^9)29 + 25^9 -1 = 0, and there is one parti- 



2 




FIG. 1: (Color online) The deformation of the crystal unit 
cell with lattice vectors ei as a function of packing fraction <j) 
in a system of hard superballs with q — 2.5. The dashed lines 
in the figures indicate the values for the Ci crystal. 



cle in the unit cell For instance, when q — 2.5, one 
finds that (61,62) = 0.5, (61,63) = (63,62) = 0.60552, 
e2|/|ei| = 1 and |63|/|6i| = 0.825737, where (6^,6^) 
is the cosine of angle between 6^ and 6j. The calcu- 
lated angles and the length ratios between lattice vec- 
tors as a function of packing fraction </) for the cube- 
like particles with q = 2.5 are shown in Fig. [T] We 
find that at packing fractions approaching close packing, 
the crystal remains in the Ci phase. With decreasing 
packing fraction, the crystal lattice deforms towards an 
fee structure: (01,62) — (61,63) = (62,63) = 0.5 and 
1^21 / |ei| = I63I / |6i| = 1. Moreover, when 1 < q < 3, it 
is found that the deformed Co and deformed Ci crystal 
melt into an fee plastic crystal phase. By Einstein in- 
tegration, we calculated the Helmholtz free energy as a 
function of packing fraction for both the fee plastic crys- 
tal and the deformed Ci/Cq crystal phases [2lj. Com- 
bined with the free-energy calculations for the fluid phase 
done by Widom's particle insertion method, we obtain 
the phase boundaries in the phase diagram shown in 
Fig. [2j The part of the phase diagram for hard cube-like 
superballs roughly a gree s with the empirical phase dia- 



gram by Batten et al. At high packing fractions, there 
are stable deformed Cq and Ci phases. When q > 1.1509, 
the close-packed structure is the Ci crystaljwhereas it 
is the Co crystal whenever 1 < q < 1.1509 [10|. To de- 
termine the location of the transition from the deformed 
Co crystal to the deformed Ci crystal, we performed two 
series of NPT MC simulations with increasing value of q 
for the first series and decreasing q for the second series 
of simulations at pressure P* = Pv/ksT ~ 250, with 
ks the Boltzmann constant, T the temperature, and v 



the volume of the particle [12[. The first series started 



from a Co crystal phase, while the second series of simu- 



lations started from a Ci crystal phase. Our simulations 
show that the phase transition occurred around q = 1.09 
at packing fraction (j) — 0.736 as shown by the asterix 
in Fig. [2] Moreover, for hard cubes {q = 00) the Ci 
crystal is a simple cubic (sc) crystal. Although it was 
found that for hard cubes there is a significant amount 
of vacancies in the simple cubic crystal, it only shifts the 
phase boundary by ^ 2% in packing fraction [l^]- In 
our simulations, we did not observe any vacancies in the 
crystals of hard superballs with g < 3, we therefore as- 
sume that the possible presence of vacancies would not 
shift the phase boundary significantly. 

The other part of the phase diagram concerns the 
octahedron-like superballs. For 0.79248 < q < 1, we ob- 
tained a denser structure than the predicted Oq lattice 
of Ref. fld\. For instance, after compressing the system 
to pressures around P* = 10'' at g = 0.85, we obtained a 
body-centered-tetragonal (bet) crystal with (f> — 0.7661. 
This is denser than the Oq crystal, which achieves cj) — 

0. 7656 a,t q = 0.85. Note however that these two crystals 
are very similar to each other, since Co is also a form of 
a bet lattice. The only difference is that the orientation 
of the particles in the Oq crystal is the same as the sym- 
metry of the axes in the crystal lattice, while in our bet 
crystal there is a small angle between these two orienta- 
tions in the square plane of the crystal. Furthermore, for 
q < 0.79248, we also found a crystal with denser packing 
than the predicted Oi crystal in Refs. [l^l- For q = 0.7, 
we performed fioppy-box MC simulations with several 
particles to compress the system to a high pressure state, 

1. e., P* = 10^. We found a deformed bee (dbcc) crystal 
shown in Fig. [21 which is an intermediate form between 
the bee lattice and the Minkowski crystal j22|. The lat- 
tice vectors are 61 = 0.912909i -I- 0.912403j - 0.912165k, 
62 = -0.271668i -I- 1.80916j - 0.288051k and 63 = 
0.28834i - 0.272001j - 1.80882k, where i, j and k are 
the unit vectors along the axes of the particle. Our dbcc 
crystal is very close the predicted Oi crystal, whose lat- 
tice vectors are ei = 0.912492i + 0.912492j - 0.912492k, 
62 = -0.2884i + 1.80629j - 0.2884k, and 63 = 0.2884i - 
0.2884j — 1.80629k. However, it has a packing frac- 
tion oi (p = 0.832839 which is denser than the predicted 
Oi crystal with = 0.824976 in Refs. [l3| by roughly 
1%. In Ref. [iq|, the Oq and Oi phases are found to 
switch a.t q — 0.79248. We also observed that the bet 
and dbcc crystals both transform into the bee phase at 
q = 0.79248. 

As shown in Fig. [21 when the shape of the superballs 
is close to spherical, i.e., 0.7 < q < 3, there is always a 
stable fee plastic crystal phase. Surprisingly, when the 
shape of superballs is octahedron-like, we find a stable 
bee plastic crystal phase. Moreover, around q — 0.8 we 
even find a fairly broad two-phase regime where a low- 
density fee plastic crystal coexists with a high-density 
bcc plastic crystal phase. In order to quantify the orien- 
tational order in the bcc plastic crystal, we calculate the 
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FIG. 2: (Color online) Phase diagram for hard superballs in the </!> (packing fraction) versus 1/q representation where q is the 
deformation parameter. Here the Ci and Co crystals are defined in Ref. [10], where the particles of the same color are in the 
same layer of stacking. The solid diamonds indicate the close packing, and the locations of triple points are determined by 
extrapolation as shown by the dashed lines. The phase boundaries for hard cubes are taken from Ref. [Q. 



cubatic order parameter 5*4 given by 11 1 
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Vf35|u, 



30 1 u. 
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where Uij is the unit vector of the j-th axis of particle i, 
N the number of particles, and n is a unit vector. The 
cubatic order parameter Sa is shown in Fig. [3] as a func- 
tion of packing fraction for a system of superballs with 
q ~ 0.7. We observe that {S^) ~ 0.2 at low packing 
fractions, which means that there is a very weak orienta- 
tional order in the system [2^. With increasing packing 
fraction, the cubatic order parameter increases monoton- 
ically to around 0.65 at a packing fraction of 0.7, which 
is indicative of a medium-ranged orientationally ordered 
system. This suggests that the entropic repulsion due to 



the rotation of the octahedron-like superballs stablizes 
the bcc lattice. 

Due to the numerical instability in the overlap algo- 
rithm, we are not able to investigate systems of super- 
balls with q < 0.7 [23f. However, we can use the sepa- 
rating axis theorem [6| to simulate hard superballs with 
q = 0.5, i.e., perfect octahedra. When we compressed the 
system from a fluid phase, we did not observe the spon- 
taneous formation of a crystal phase in our simulation 
box within our simulation time. When we expand the 
Minkowski crystal, which is the close-packed structure 
of octahedra, in NPT MC simulations by decreasing the 
pressure, the system melts into a bcc plastic crystal phase 
as shown in Fig. S) We also calculated the free energy 
for these three phases to determine the phase boundaries. 
To exclude finite-size effects in the free-energy calcula- 
tion of crystal phases, we performed Einstein integration 
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FIG. 3: (Color online) Cubatic order parameter S4 as a func- 
tion of packing fraction (f> for a bee plastic crystal phase of 
hard superballs with q — 0.7. The inset shows a typical config- 
uration of a bcc plastic crystal of hard superballs with q = 0.7 
at ^ = 0.54. 

for systems of = 1024, 1458, and 2000 particles, and 
applied a finite-size correction [2l|. We confirmed the 
errors in the free-energy calculations to be on the or- 
der of 10~^A;Br per particle. The calculated free-energy 
densities for the three phases are shown in Fig. S) Em- 
ploying a common tangent construction, we found that 
there is only phase coexistence between a fluid phase and 
a Minkowski crystal phase, while the bcc plastic crystal 
phase is metastable. However, the free-energy differences 
between the fluid and the plastic crystal phase at the bulk 
coexistence pressure is very small , i.e., ~ 10~^kBT per 
particle. Moreover, the Minkowski crystal melts into a 
bcc plastic crystal before melting into the fluid phase. 
Our results thus show that the rounded corners of octa- 
hedra play an important role in stablizing the bcc plastic 
crystal phase which is a new plastic crystal phase for 
systems of hard particles. 

In conclusion, using free-energy calculations we have 
determined the full phase diagram of hard superballs 
with shapes interpolating between cubes and octahedra, 
i.e., 0.5 < q < 00. In systems of cube-like superballs 
{q > 1), we find a stable deformed Ci phase at high 
packing fraction, except close to the sphere- limit {q = 1) 
where a deformed Co crystal is stable. For q < 3 the crys- 
tal phase melts into an fee plastic crystal before melting 
into a fluid phase of cubelike superballs. In systems of 
octahedron-like superballs (0.5 < g < 1), we flnd a stable 
bet or a deformed bcc crystal phase upon approaching 
close packing, with a crossover at q = 0.79248. More- 
over, a stable fee plastic crystal appears at intermediate 
densities for 0.7 < g < 1. Interestingly, for q < 0.85, we 
find a novel stable bcc plastic crystal phase, which can 
even coexist with the fee plastic crystal phase at around 



FIG. 4: (Color online) A part of the equation of state for 
hard octahedra. The pressure Pv/ksT and free-energy den- 
sity F/V — pfjLc -I- Pc as a function of packing fraction (j>. Here 
V is the volume of the particle; F and V are the Helmholtz 
free energy and the volume of the system (in units of parti- 
cle volume) respectively; /ic and Pc are the chemical potential 
and pressure at bulk coexistence respectively with p the num- 
ber density of the particles. The solid lines in the EOS for 
the Minkowski and the bcc plastic crystal phases are obtained 
by melting the close-packed Minkowski crystal in floppy box 
NPT MC simulations, and the dotted line for the bcc plas- 
tic crystal is obtained by compressing the crystal in cubic box 
NPT MC simulations. The black points and dashed line show 
the coexistence between the fluid phase and the Minkowski 
crystal phase. 



q — 0.8. More surprisingly, the bcc and fee rotator phases 
are unstable for the flat-faced and sharp-edged hard octa- 
hedra and hard cubes, respectively, which suggests that 
the rounded corners of superballs may play an important 
role in stabilizing rotator phases. It is tempting to argue 
that entropic directional forces that tend to align 
sufficiently large flat faces of polyhedral-shaped particles 
destabilize rotator phases in favor of crystals. We stress 
here that rounded corners are not a necessary condition 
for stable rotator phases since almost spherical polyhe- 
dral particles have been shown to form rotator phases 
as well Finally, we also observed a two-step melting 
phenomenon in the system of hard octahedra, such that 
the Minkowski crystal melts into a metastable bcc plastic 
crystal before melting into the fluid phase. Nanoparticle 
self-assembly is therefore surprisingly sensitive to particle 
curvature. 
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and Utrecht University High Potential Programme. 



5 



* Electronic address: r.ni@uu.i;l| 
^ Electronic address: m.dijkstral@uu.nl| 
[1] S. C. Glotzer and M. J. Solomon, Nature Mater. 6, 557 
(2007). 

[2] A. Yethiraj and A. van Blaaderen, Nature (London) 421, 
513 (2003). 

[3] Y. Zhang, F. Lu, D. van der Lelie, and O. Gang, Phys. 

Rev. Lett. 107, 135701 (2011). 
[4] A.-P. Hynninen, J. H. J. Thijssen, E. C. M. Vermolen, 

M. Dijkstra, and A. van Blaaderen, Nature Mater. 6, 202 

(2007). 

[5] A. Haji-Akbari, M. Engel, A. S. Keys, X. Zheng, R. G. 

Petschek, P. Palffy-Muhoray, and S. C. Glotzer, Nature 

(London) 462, 773 (2009). 
[6] U. Agarwal and F. Escobedo, Nature Mater. 10, 230 

(2011). 

[7] S. Torquato and Y. Jiao, Nature (London) 460, 876 
(2009). 

[8] T. C. Hales and S. P. Ferguson, Discrete Comput. Geom. 
36, 5 (2006). 

[9] J. de Graaf, R. van Roij, and M. Dijkstra, Phys. Rev. 

Lett. 107, 155501 (2011). 
[10] Y. Jiao, F. H. Stillinger, and S. Torquato, Phys. Rev. E 

79, 041309 (2009); ibid, Phys. Rev. E 84, 069902 (2011). 
[11] R. D. Batten, F. H. Stillinger, and S. Torquato, Phys. 

Rev. E 81, 061105 (2010). 



[12] M. Marechal, A. Cuetos, B. Martmez-Haya, and M. Di- 
jkstra, J. Chem. Phys. 134, 094501 (2011). 

[13] F. Smallenburg, L. Filion, M. Marechal, and M. Dijkstra, 
submitted (2011), arXiv:1111.3466v2 [cond-mat.soft]. 

[14] C. Vega and P. A. Monson, J. Chem. Phys. 107, 2696 
(1997). 

[15] M. Marechal and M. Dijkstra, Phys. Rev. E 77, 061405 
(2008). 

[16] P. Bolhuis and D. Frenkel, J. Chem. Phys. 106, 666 
(1997). 

[17] P. G. Bolhuis, D. Frenkel, S.-C. Mau, and D. A. Huse, 
Nature 388, 235 (1997). 

[18] L. Rossi, S. Sacanna, W. T. M. Irvine, P. M. Chaikin, 
D. J. Pine, and A. P. Philipse, Soft Matter 7, 4139 (2011). 

[19] L. Filion, M. Marechal, B. van Oorschot, D. Pelt, 
F. Smallenburg, and M. Dijkstra, Phys. Rev. Lett. 103, 
188302 (2009). _ 

[20] See Supplemental Material at http: //unknown for tech- 
nical details and the movie of plastic bcc crystal formed 
by superballs with q — 0.7 at packing fraction 0.54. 

[21] D. Frenkel and B. Smit, Understanding Molecular Sim- 
ulation: From Algorithms to Applications (Academic 
Press, 2002). 

[22] H. Minkowski, Konigliche Gesellschaft der Wis- 
senschaften zu Gottingen pp. 311-355 (1904). 

[23] A. Donev, Ph.D. thesis, Princeton University (2006). 

[24] P. F. Damasceno, M. Engel, and S. C. Glotzer (2011), 
arXiv: 1 1 09 . 1323v2 [cond-mat . soft] . 



